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Abstract 

In this paper, we calculate the branching ratios and the direct CP- violating asymmetries for 
decays B° -» ag(980)vr , a£(980)7i-, (980)vr+ and B~ -» ag(980)vr-, (980)vr° by employing 
the perturbative QCD factorization approach. Although the light scalar meson ao(980) is widely 
perceived as primarily the four-quark bound states, the calculation shows that the 2-quark model 
supposition for ao(980) can not be ruled out by the currently available experiment upper limits. 
In these considered decays, the branching ratio for the B~ — > aQ(980)7r~ is the largest, 2.8 x 10 -6 , 
while its direct CP asymmetry is the smallest, ~ 14%. Although the direct CP asymmetries for 
the decays B° -> a[](980)7r°, B~ -> 00(980)^° are large, about (70 ~ 80)%, it is still difficult to 
measure them, since their branching ratios are small, around (4 ~ 5) x 10~ 7 . 
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I. INTRODUCTION 



The study about scalar meson is an interesting topic for both theory and experiment. 
In order to uncover their mysterious structure, intensive studies have been done for the 
B meson decays involving a scalar meson as one of the two final state mesons. Such kind 
of decays have been studied by employing various factorization approaches, such as the 
generalized factorization approach the QCD factorization (QCDF) approach 0, Hj|, 
the perturbative QCD (pQCD) approach 0,1,0,0,9, and by using the QCD sum rule 

1 03.01. 

On the experimental side, from the first scalar meson /o(980) observed in the decay 
mode B — > f (980)K by Belle 12J, and confirmed by BaBar [13( later, then many channels 
involving a scalar in the final state have been measured by Belle 14], ll5| and BaBar 



16, [l7j. For example, the decays B — > ao(980)7r were searched by BarBar five years ago 



181 ]. especially the decay B — > a (980)7r°, which has been considered as a best candidate 



to distinguish the nature of the scalar a (980) [19J. The authors of Ref.jl9j argued that 
if the branching ratio of this channel can be measured accurately by the experiment, 
one can separate the four- and two-quark assignments, because the results of these two 
assignments have a difference of one order of magnitude. So in the past three years, 
BarBar have given this channel twice measurements [lo|, 2l[ and get two almost identical 
upper limits. For our considered decays, only the experimental upper limits are available 
now for some of them 22]: 



Br{B° 
Br(B- 
Br(B~ 
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c4(980)tt-) < 3.1 x 10 

• ag (980)tt°) < 1.4 x 1(T 6 

* ag(980)7r _ ) < 5.8 x KT 6 



In this paper, we will study the branching ratios and CP asymmetries of B° — > 
ag(980)7r°,4(980)7rT and B~ -> a (980)vr°, ag(980)7r~ within perturbative QCD ap- 
proach based on kr factorization. In the following, we use ao to denote ao(980) in some 
places for convenience. It is organized as follows. In Sect 0U the status of the study on 
the physical properties of ao, the relevant decay constants and light-cone distribution am- 
plitudes are discussed. In SecJTTIl we then analysis these decay channels using the pQCD 
approach. The numerical results and the discussions are given in the section IIVI The 
conclusions are presented in the final part. 



II. PHYSICAL PROPERTIES OF THE FINAL PARTICLES 

Many scalar mesons below 2GeV have been found in experiments. We cann't ac- 
commodate these scalar mensons into one nonet, but need at least two nonets be- 
low and above 1 GeV[23T]. Among them, the scalar mesons below 1 GeV, including 
/o(600)(a),/ (980), 1^(800) (/c) and a (980), are usually viewed to form an SU(3) nonet; 
while scalar mseons above 1 GeV, including / (1370), / (1500)// (1700), if*(1430) and 
ao (1450) form the other SU(3) nonet. There are several different scenarios to describe 
these mesons in the quark model 0, 25. l26l . [27|. For example ao(980) meson, which has 



been suggested as qq lowest lying state [24J (called scenario I) or four-quark bound state 
[25[ (called scenario II). In the scenario I, the former SU(3) nonet mesons are treated as 
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the qq ground stats, and the latter nonet ones are the first excited states; in the scenario 
II, the former nonet mesons are viewed as four-quark bound states, while the latter nonet 
ones are qq ground states. Some people also consider that it is not made of one simple com- 



ponent but might have a more complex nature such as having a KK component [26J, 127] . 
even the superpositions of the two- and four- quark states. In order to make quantitative 
predictions, we identify ao(980) as the two-quark state in the calculation. 
In 2-quark model, the decay constants for scalar meson ao are defined by: 

(oo(p)|?27m9iI°) = /ao2V> < a o(p)|^2<?i|0> = m a J ao . (2) 

For the neutral scalar meson ao cannot be produced via the vector current (restricted by 
the charge conjugation invariance or the G parity conservation), the vector decay constant 
f ao = 0. As to the charged scalar mesons ag , from the equation of motion: 

m a- 

U„-f n - = /„-, with a - = — — q — , (3) 

its vector decay constant is proportional to the mass difference between the constituent 
u and d quarks. It is easy to see the vector decay constant is very small, and will equal 
to zero in the SU(3) limit. So we only need consider the scalar decay constant f ao , which 
is scale dependent. Fixing the scale at 1 GeV, the value is f ao = (365 ± 20)MeV, which 
is calculated in QCD sum rules 

The light-cone distribution amplitudes (LCDAs) for the scalar meson ao can be written 

as: 

(a (p)|?i(z) i g 2 (0) i |0) = -±= [ dxe^ 

\?2,l\ c Jo 

x fcteao 0*0 + m a $U X ) + m *o ( i+i- - 1) $r o ( X ) }jl , ( 4 ) 

where n + and n_ are light-like vectors: n + = (1, 0, Oy), n_ = (0, 1, Ot), and n+ is parallel 
with the moving direction of the scalar meson ao- The normalization can be related to 
the decay constants: 



J « ao (x) = J dx$l(x) =0, J dx$ s ao (x) = 
The twist-2 LCDA can be expanded in the Gegenbauer polynomials: 



(5) 



$ ao (x,/i) = — =/ ao (/i)6x(l -x)J2 B m (v)C 3 J 2 (2x - 1), (6) 

" c m=l 

the values for Gegenbauer moments Bi, B 3 have been calculated in [^] as: 

B v = -0.93 ±0.10, B 3 = 0.14 ±0.08. (7) 

These vaues are taken at jj, — 1 GeV and the even Gegenbauer moments vanish. 

As for the twist-3 distribution amplitudes <3>^ o and $J o , they have not been studied in 
the literature, so we adopt the asymptotic form 



1 /• ,!,/ 1 



$ao 2v^Ar /ao ' $ao 2v/2iv/ ao(1 2X) ' 
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III. THE PERTURB ATIVE QCD CALCULATION 



Under the two-quark model for the scalar meson a supposition, we would like to 
use pQCD approach to study B decays into a and n. The decay amplitude can be 
conceptually written as the convolution, 

A(B -> 7ra ) ~ Jd%d%d% Tr [C^)M*i)M*2)$«o k 2 , k 3 , t)] , (9) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C (t) is the Wilson coefficient which results from the radiative 
corrections at short distance. In the above convolution, C(t) includes the harder dynamics 
at larger scale than Mb scale and describes the evolution of local 4-Fermi operators from 
raw (the W boson mass) down to t ~ 0(^/AMb) scale, where A = M#— m^. The function 
H(ki, k 2 , kz, t) describes the four quark operator and the spectator quark connected by a 
hard gluon whose q 2 is in the order of AMg, and includes the 0(\/AMb) hard dynamics. 
Therefore, this hard part H can be perturbatively calculated. The function $(7r,a ) are 
the wave functions of tt and do, respectively. 

Since the b quark is rather heavy we consider the B meson at rest for simplicity. It is 
convenient to use light-cone coordinate (p + ,P~,Pt) to describe the meson's momenta, 

p ± = ^=(p°±p 3 ), and Pt = (p\p 2 )- (10) 

Using these coordinates the B meson and the two final state meson momenta can be 
written as 

P B = ^|(l,l,O r ), P 2 = ^|(1,0,0 T ), P 3 = ^|(0,1,0 T ), (11) 

respectively. The meson masses have been neglected. Putting the anti- quark momenta 
in B, P and S mesons as k%, k 2 , and k 3 , respectively, we can choose 

ki = (xxP+^Mt), k 2 = (x 2 P+,0,k 2T ), k 3 = (0,x 3 P 3 -,k3 T ). (12) 

For these considered decay channels, the integration over A;i 5 ; ciiici kn in eq.([n]) will 
lead to 



A(B -> 7ra ) ~ ji 



dxidx 2 dx ?J bidbib 2 db 2 b 3 db 3 

■Tr [C(t)$ B (x 1 ,6 1 )$ w (x 2 ,6 2 )$ ao (x3,63)^(^,^,t)5t(^)e- s W] , (13) 

where bi is the conjugate space coordinate of kix, and t is the largest energy scale in 
function H(xi,bi,t). In order to smear the end-point singularity on Xj, the jet function 
S t {x) [28], which comes from the resummation of the double logarithms ln 2 Xj, is used 
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where the parameter c = 0.4. The last term e~ s ^ in Eg. (1131) is the Sudakov form factor 
which suppresses the soft dynamics effectively (29| . 

For the considered decays, the related weak effective Hamiltonian i? e // can be written 

as 



10 

(dGu)0?(//) + C 2 (//)OfGu)) diri Oi(M) 



H ef f = -j= V lb V qd 



V2 

q=u,c 



=3 



(15) 



with the Fermi constant Gf = 1.16639 x 10 5 GeV 2 , and the CKM matrix elements V. 
We specify below the operators in H. e ff for b — > d transition: 

= da^Lu/3 ■ upypLbo, , 0% = d_ a ^Lu a ■ u^^Lbp , 

3 = d_ a y»Lb a ■ J2 q > q'pl» L q'p , °i = da^Lbfs ■ J2 q > Qpl^ L Q'a > 

5 = d a YLb a ■ J2 q < Q.'pl^Ro'p , 6 = d a y>*Lbp ■ q' p l^Rq' a , (16) 

7 = Idt/fLba ■ i2 q > e q'q.'pln R( l'p , °s = Ida^Lbp ■ Y^ q , e ql q'^^Rq' a , 

9 = ld Q YLb a ■ J2 q > e q >q'^»Lq'i3 , °io = \d a ^Lb p ■ J2 q > ^Qpl^Wa > 

where a and j3 are the SU(3) color indices; L and R are the left- and right-handed 
projection operators with L = (1 — 75), R = (1 + 75). The sum over q' runs over the quark 
fields that are active at the scale \i = 0(mb), i.e., (q'e{u, d, s, c, b}). 

In the following, we take the B° — > 7r°aQ decay channel as an example to expound. 
There are 8 type diagrams contributing to this decay, as illustrated in FigJTJ For the 
factorizable emission diagrams (a) and (b), operators Oi_4 )9) i are (V— A)(V — A) currents, 
and the operators Os_g have a structure of (V — A)(V + A) , the sum of the their amplitudes 
are written as F en and F^ 1 . In some other cases, we need to do Fierz transformation for 
the (V — A)(V + A) operators and get (S — P)(S + P) ones which hold right flavor 
and color structure for factorization to work. The contribution from the operator (5 — 
P)(S + P) type is written as F^ 2 . Similarly, for the factorizable annihilation diagrams 
(g) and (h), the contributions from (V - A)(V - A),(V - A)(V + A),(S - P)(S + P) 
currents are F a7T ,F^ and F^ 2 . For the nonfactorizable spectator diagrams (c, d) and 
the nonfactorizable annihilation diagrams (e, f), these three kinds of contributions can 
be written as M en , M^, M^ 2 and M aw , M^J^M? 2 , respectively. Since these amplitudes 
are similar to those B -> f (980)K(-K, 77W) @, [J or B -> 00(980)^ [8], we just need to 
replace some corresponding wave functions and parameters. 

Combining the contributions from different diagrams, the total decay amplitudes for 
these decays can be written as: 

2M(B° -> a°vr°) = [(-M e7r + M an + M eao + M aao )C 2 + (F aw + F eao + F aao )a 2 ] 

-6 |W + m£ + m£ + MZ)(C 5 - ic r ) 

+ (M en + M an + M eao + M aao ) (C 3 + 2C 4 - + ~C 10 ) 

+ (Mi; 2 + Mf w 2 + Mf a 2 + Mf a 2 )(2C 6 + l -C,) + (F a7r + F eao 

111, 
+Faao) (2a 3 + a 4 - 2a 5 - -a 7 + -a 9 - -a w ) 

+ K 2 + F P J + F e P2 + F a P2 ) (a 6 - las) } , (17) 
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FIG. 1: Diagrams contributing to the decay B° — > vr°aQ . 



_M(5 -> a vr H 



£ u [^2 + M e7r C a + M aao C 2 ] - 6 j^C \C S - -C 7 
+M e p n \C 5 + C 7 ) + M a7r (C 3 + C 4 - ic 9 - ic 10 ) + M aao (C 4 + C 10 ) 

M e7r (C s + C») + M a C 2 (C 6 - ^C 8 ) + M2(C 6 + C 8 ) 
111 

+-fa7r(«3 + 0-4 — + -07 — -Og — ^W) 



+F aao (a 3 + a 9 - a 5 - o 7 ) + F e7r (a 6 + a 8 ) + F a7r ( a 6 - -o 8 ) J> , (1£ 



^(fi -> a+O = [i^eaoai + F an a 2 + M eao d + M aw C 2 ] - & |m£ ( C 5 - ~C 7 

+M£(tf 5 + C 7 ) + M aao (C 3 + C 4 - J<7 9 - ic? 10 ) + M aw (C 4 + C 10 ) 



2 2 

rP2, 



M eao (C 3 + C 9 ) + M£ (C 6 - -C 8 ) + M a ^(C 6 + 0.) 
+F aao (a 3 + a 4 — a 5 + -07 — -a 9 — t^io) + -^ea (a-4 + a-10) 
+F an (a 3 + a 9 - a 5 - a 7 ) + F P2 (a 6 + a 8 ) + ( a 6 - -a 8 ) L (19) 



TABLE I: Input parameters used in the numerical calculation 



Masses 


m ao = 


0.9847 GeV, 




= 1.3 GeV, 




M B = 


= 5.28 GeV, 




= 0.14 GeV, 


Decay constants 


fB = 


0.19 GeV, 


U 


= 0.13 GeV, 


Lifetimes 


T B ± = 1. 


671 x 10~ 12 s, 




1.530 x 10" 12 s, 


CKM 


Vtb 


= 0.9997, 


v td = 


0.0081e- i21 - 6 °, 




Vud 


= 0.974, 


v ub = 


0.00393e- i60 ° . 



V2M(B- -> a° n-) = £ u [M en C 2 + (-M a7r + M eao + M aao )C t + {-F a7T + F eao + F aao ) ai ] 

-6 {-M£{C B - \c 7 ) + (-ilC 1 + M£ + M£)(C 5 + C 7 ) 



2 

■ ^ I, n ! + \C 9 + ^C 10 ) + (-M a7r + M eao + M aoo ) (Cg + Cg) 



3 1 

+ -C 8 M£ 2 + (-F a7r + F eao + F aao ) (a 4 + aio) - F™(a 6 - -a 8 ) 

+ (-^? + ^ao+^o)K + a 8 )}, " (20) 



V2M(B~ -> a tt°) = £ u [M eao C 2 + (-M aao + M e7r + M a7r )d + F e a a 2 + ( — 

1 

— I 

2 



-6 <! -aC(c 5 - -c 7 ) + (-mZ + m p J + m£)(c s + c 7 ) 



+M eao (-C 3 + -Cg + ^Cio) + (-M aao + M e7T + M a7T ) (C 3 + C 9 ) 



+ -C 8 Mf a 2 o + (-F aao + F en + F a7T ) (a 4 + a 10 ) - F P2 (a 6 - -a 8 ) 
+F eao (-a 4 - ^a 7 + + ±a 10 ) + (-F a p a 2 + i^ 2 + F P2 )(a 6 + a 8 ) |> , (21) 



where £ u = V u bV* d , £ t = V t bV t * d . The combinations of the Wilson coefficients are defined 
as usual |3l|: 



ai(/i) 


= cm + 


cm 

3 ' 




= cm + 


3 




= aw + 


3 



z — 3 ; 5 ; 7 ; 9 ; 



1 = 4,6,8,10. (22) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In the numerical calculation, we will use the input parameters as listed in Table [I] 
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TABLE II: Branching ratios (xlCT 6 ) for the decays B° -> agvr , a^vr^ and B~ -> a Q vr°, a%n~. 
The first theoretical error is from the the scalar meson decay constant, the second and the third 
one are Gengebauer moments B\ and B3 for twist-2 LCDAs of a(980). 



Channel 


This work 


Data 


QCDF [3] 


B° -» a° 7r° 


n r-i +0.08+0.09+0.00 
u -' Ji -0.07-0.09-0.00 




0.2 


B° -» a^ir- 

B° -» CLQ7T+ 
B° -> Oq 7T~ + ag 7T+ 


n Qc+o. 10+0. 14+0.01 

u -°°-0.09-0.14-0.00 
n c-i +0.05+0.09+0.07 
u ' di -0.06-0.09-0.06 
n Qo+0.10+0.15+0.02 
u - JJ -0.10-0.14-0.00 


< 3.1 


7.6 
0.6 


B~ ao7r° 
B- -» agvr" 


n A i+o.oo+o.oo+o.oo 

u -^ i -0.13-0.14-0.12 
9 Q+0.00+0.00+0.00 
z -°-0.79-0.85-0.58 


< 1.4 

< 5.8 


0.2 
3.4 



In the B-rest frame, the decay rates of B — > a (980)7r can be written as: 

r =3Jy^ |2(i - ry (23) 

where r ao = m ao /mB and .M is the total decay amplitude of B — > a (980)7r, which has 
been given in section IIHI 

Using the wave functions as specified in previous section and the input parameters 
listed in Table [U it is straightforward to calculate the CP-averaged branching ratios for 
the considered decays, which are listed in Table [TT1 In this table, we have included 
theoretical errors arising from the uncertainties in the scalar meson decay constant f ao 
and the Gengebauaer moments B\ and B% for twist-2 LCDAs of ao(980). 

From the numerical results, one can find that: Firstly, the branching ratio of B° — > a^ 
is larger than that of B° — > /o(980)7r° [?}, for the small uu and dd component in the 
/o(980), but much smaller than the branching ratio Br(B° — > 7r°7r°) = (1.62±0.31) x 1CT 6 . 
Since the scalar meson ao(980) has vanishing decay constants in the isospin limit, the ao~7r + 
rate is smaller than that of its conjugated channel. Secondly, the ratio of these two rates is 
about 1/2. If one includes the interference between these two decay modes, the branching 
ratio of B° — > a '7r + + a^vr - is close to 1 x 10~ 6 . Thirdly, for the two charged decays, it 
is the same reason with the previous two neutral ones, and the aQ~7r rate is much smaller 
than that of clqTT^ ■ Lastly, the 2-quark model supposition of ao(980) can not be ruled out 
by the current experimental data, and this point is different from the QCDF prediction 
[3|. Since only the upper limits for these channels are available now and the daughter 
branching fraction has been taken to be 100%, so the further accurate data are needed to 
clarify this discrepancy. 

In Fig. [21 we plot the branching ratios of B° — > a l 7r°,aQ7r + + a^-K^ and B~ — > 
aQ7r~,ag7r° as functions of the CKM angle a = arg[— v i d y b J . From these figures, it is 

found that the branching ratio of B° — > Oq tt + + a^ir^ is more sensitive to the angle a 
than those of other three decay channels. It has been claimed in Ref. 0] that the B° — B° 
interference plays no role in the a^n T channels. We disagree and we argue that though 
the B° — B° interference in the a^ir T channels may be not as strong as that in some other 
channels, such as ai(1260) =F 7r ± modes. There must exist certain interference among them, 
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FIG. 2: Branching ratios (in units of 10 -6 ) of (a) B° — ► a^tr (solid curve), B° — > ag 7r + + aQ~7r _ 
(dashed curve) and (b) £>~ — > <2q 7r° (solid curve), £?~ — > a^-K~ (dashed curve) decays as a 
function of the CKM angle a. 



as can be seen easily in the numerical values in Table II, so it may be unreasonable to 
regard the decays B° — > a^ii^ as the self-tagging decays in the experiment. 

Now we turn to the evaluations of the CP-violating asymmetries of B — > Oq 7r°, o^ti" 
and B a^it , a^n T decays in pQCD approach. For the charged decay channels, the 
direct CP-violating asymmetry can be defined as: 

For the neutral decays B° — > a^ir , there are both direct CP asymmetry Afj P and 
mixing- induced CP asymmetry A™p . The time dependent CP asymmetry of B decay 
into a CP eigenstate / (aQ7r ) is defined as: 

Acp(t) = A d £ P {B d -> /) cos(AMt) + Atp{B d -> /) sin(AMt), (25) 

with 

•AcpiBd /) = , -4g? (B„ -> /) = , (26) 

where 77 = ±1 depends on the CP eigenvalue of /, AM is the mass difference of the two 
neutral B meson eigenstates. 

As to the decays B — > a^ir T , since both B° and B° can decay into the final state 
(1q7i~ and the four time-dependent decay widths for B°(t) — > ao7r~,£>°(t) — > 

Oq 7r + , £>°(t) — > ag 7r + and — > Oq 7r~ can be expressed by four basic matrix elements: 

5 = (a+7r|# e// |£°), /i = (a+7r-|# e// |S°>, 

5 = (ao K + \H eff \B°), h = (a n + \H eff \B°). (28) 
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FIG. 3: The direct CP asymmetries (a) of the decays B° — > a[]vr (solid curve), B~ — > a^n 
(dashed curve), B~ — > a^ir" (dotted cureve) and the CP asymmetry parameters (b) of the 
decay B° — > a^~7r — + aQir + : a e / (dash-dotted curve) , (dotted cureve), a e+e i (dashed cureve) 
and a e+? / (solid curve) as functions of the CKM angle a. 



Following the notation of Ref. 32|, |33j, the four CP violating parameters are given by 
the following formulae: 

-2im(5|) 

le+e ' = l + \h/g 9 \ 2 
-2/m(;f) 

le+ " 1 + \9lh? ' 



a e i = ■ 


9 


2 _ 


h\ 2 


9 


2 + \h\ 2 





9? 


\h 


2 + 


9? 



where q/p = e 21/3 is denned via 



B l = p\B {, ) + q\B ) 



Bo 



p\B°)-q\B"), 



(29) 



(30) 



with \p\ 2 + \q\ 2 = 1 and (3 being a CKM angle. 



From the Figj3](b), one can find the central values of the CP- violation parameters: 

0.31, a e+e > = 0.94, a- e , = 0.93, a t+ - t > = 0.32; (31) 



for a = 100°. 

The direct CP asymmetries are shown in Figj3]^a). The branching ratio of decay B~ — > 
ag7r~ is the largest one among the considered channels, while its direct CP asymmetry is 
the smallest, 14%. Although the CP asymmetries for the decays B° — > a^(980)7i° , B^ — > 
<2q (980)7t° are large, about (70 ~ 80)%, it is still difficult to measure them, since their 
branching ratios are small, around (4 ~ 5) x 10~ 7 . 



V. CONCLUSION 

In this paper, we calculate the branching ratios and CP-violating asymmetries of B° — > 
ciqTT , B° — > (Iq7t~, B° — > ciq7t + , B~ — > (Iq7t~ and B~ — > ci^tt decays in the pQCD 
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factorization approach by identifying ao(980) as the 2-quark content. Using the decay 
constants and light-cone distribution amplitudes derived from QCD sum-rule method, we 
find that: 

• Since the scalar meson ao(980) has vanishing decay constant in the isospin limit, one 
can find Br(B° -> a+vr ) > Br(B° -> a ir + ) and Br(B~ -> a°7r~) > Br(B~ -> 
a -7r°). 

• If one includes the B° — B° interference in the a^n T channels, Br(B° — ► a[j~7r _ + 
a -7T + ) = 0.93 x 10~ 6 , which is close to £r(£° -> a^Tr") = 0.86 x 10~ 6 . There must 
exist certain interference in these two neutral channels. 

• Although the CP asymmetries of B° — > aQ(980)7r°, B~ — > 0^(980)^° are large, 
about (70 ~ 80)%, it is still difficult to measure them, since their branching ratios 
are small. 

• For our predictions agree well with the currently available experimental upper limits, 
the 2-quark model supposition for a (980) can not be ruled out. This point of view 
is different from the conclusion obtained by using the QCDF approach. 
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